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Abstract 

Quantum dynamics of integrable systems is discussed. Localized wave 
packets generalizing the conventional coherent states of minimal uncer- 
tainty are constructed. The wave packet moves along a certain trajectory 
and does not change its shape for times of order i. 

In this Letter we suggest that the definition of a coherent state should be re- 
lated to the Hamiltonian. The conventional coherent states jl| are defined as 
eigenstates of the annihilation operator 

a | z >= z\ z > (1) 

The real and the imaginary parts of the complex variable z can be expressed by 
the mean values of the position x and momentum p operators in the coherent 
state. Under the oscillator Hamiltonian evolution Ut 

U t \z >= \exp(~iuot)z > (2) 

where u> is the oscillator's frequency. Then, because the state | z > is con- 
centrated around x = ^J^Rez it follows that Ut\z > has its support on the 

classical trajectory x(t) = ^J~^-Re(exp(~iojt)z). The simple classical evolution 

z — > exp{~iujt)z is a consequence of the simple representation of ^fujx + i-^== = 

\flexp{id) in terms of the action-angle variables (/, 9) for the harmonic oscilla- 
tor. 

We are going to generalize the construction of coherent states to integrable 
non-linear systems. It is known that the quantum non-linear dynamics cannot 
be transformed into a coherent state dynamics exactly (see the discussion in 
[D||;these authors define generalized coherent states, which do not coincide 
with ours). However, wave packets approximately localized around the classical 
trajectory of an electron in the hydrogen atom have been constructed Q . 
For this purpose a large dynamical symmetry relating the hydrogen atom to the 
fourdimensional oscillator has been utilized. In this Letter we do not assume 
any symmetry, but the integrability of a classical Hamiltonian system. Our 
construction of wave packets can be applied to multidimensional anharmonic 



oscillators. The localized states can describe a classical behavior of quantum 
molecules ||. 

The wave packets are sums over energy eigenstates ^p r {x) (x G R n ; here 
and later on we avoid the vector notation and omit vector indices if there is no 
danger of confusion). For anharmonic oscillators the leading behavior for large 
|x| of ijj r (x) is determined by the behavior of the ground state X @ We write 
the ground state in the form x = ex p{—j-) ■ It is the solution of the Hamiltonian 
eigenvalue problem with the eigenvalue E q 

h 2 

H x =(- — A + V)x = E qX (3) 
Let us make the similarity transformation 

H = TT 1 x~ 1 H x = -^-& + — VSV (4) 
2m m 

Then, H has a formal limit when Tj, — > equal to ^VSV (note that S is real if 
we restrict ourselves to systems with a non-degenerate ground state x )■ 

Let V be an analytic function on the configuration space. Let us consider the 
Heisenberg equations of motion on a state exp(—j-)V ( me mark position and 
momentum operators by hats in order to distinguish them from an argument of 
a wave function) 

dxt . S _ / S , * 
- m—exp(- — )V = ptexp(- — )V (5) 

where 

x t = U t xU+ (6) 
Using eq.(4) we can rewrite eq.(5) in the form 

_ m ^l exp (- Yp = exp (- )uJ-iW + iVS)U+V (7) 
at n n 

where Ut is the unitary group generated by H (eq.(4)). The limit h — * in 
eq.(7) can be obtained explicitly. This is so because hV acting on an analytic 
function of ti gives a function of order h. Then, in the limit Ti — > 

U t V(x)=V(x t (x)) + 0(h) (8) 

where Xt (x) is the solution of the c- number equation 

TO ^± = -iVS(x t ) (9) 

with the initial condition x. 

Neglecting terms of order fi we can see that equation (7) is simplified to 

-m^^exp(-^r)V = exp(-^r)(i\7S(x t ))V 
at n n 
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Hence, by iteration we obtain the equality 

f(x t )Vexp(~) = f( Xt )Vexp(~) (10) 

for an arbitrary analytic function f, where x t (without the hat) is the solution 
of the c- number differential equation (9). 

The definition of the coherent state \z > is now determined by the require- 
ments 

i) |z = 0>= X 

ii) \z > is transformed into a certain \z(t) > under the quantum Hamiltonian 
evolution (8) in the limit h — * 0. 

For the construction of \z > let us consider a solution S^(x) of the "imagi- 
nary" Hamilton- Jacobi equation (outside the classical domain) 

-^(VS I cl f + V = E cl (I) (11) 

depending on n parameters Ik which can be related to the action variables. We 
assume that E c i = E c i(I = 0) is the classical ground state energy. It is known 
that E q — E c i ~ 0(h). We denote by S c i the solution of eq.(ll) with the ground 
state energy E c \. The difference between S and S c i is negligible in the limit 
h — > 0. This can be seen from equation (3) which when expressed by S reads 



^Sf + ^AS + V = E q (12) 



Subtracting eqs.(ll) and (12) we can conclude that S — S c i — 0(h). 

The search for a state \z > fulfilling our assumptions i)-ii) is simplified if 
we find operators Qk which under the quantum evolution exp(-itH) transform 
multiplicatively 

Qk\ip>^ \k(t)Q k \ip> (13) 
on a dense set of vectors \ij)>. We find that the proper guess is 

Qk - exp(^\ I=0 (14) 

From eq.(8) and the estimate on S — S c i it follows that H in the limit % — > 
generates the complex classical evolution 

ft - -s v *"<« (15) 

We show that in this limit the variables Qk(t) fulfil linear equations. In fact, 
from eq.(14) 

dQk (c s —i d 2 S c i dS c i 

—rriZt) = — or M -rr^Qkdt) (16 

dt m OlkO^r O^r 
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Then, from eq.(ll) it follows that 

1 d 2 S cl dS cl dE c , 



m dx r dl s dx r dL 
Hence, setting / = in eq.(17) 

dQ k 



(17) 



dt 

where 



(0 = -iukQk® (18) 



u* = ^(/ = 0) (19) 

It follows that in eq.(13) A&(£) — exp(—iuj k t). 

Now, in terms of the variables Q k our definition of a coherent wave packet 
reads 

< x | z >= exp{~ — + W - ^ y/uJkZkQk) (20) 

fe=i 

The linear dependence on Q of the exponential in eq.(20) has as a consequence 
that \z > is an eigenfunction of a generalized annihilation operator 



,Ji d_ 1 g5 



a fc 2 >= z k \z > 



Namely 
It follows that 

< z I ^("fe +a+)\z><z\z >" 1 = i?ez fc (22) 

and 

< z \ —(a k ~ a+) \ z >< z \ z >" 1 = Jmz fe (23) 

In this way we obtain an interpretation of z as a mean value of a certain 
observable. Furthermore, the eigenvalue condition is a necessary condition for 
a minimal uncertainty property of the operators which constitute the real and 
imaginary part of a (see refs.|| and jHj). It follows from eqs.(8),(18) and (20) 
that under the quantum evolution 

U t \z >= \exp(~iu;t)z > +0(H) 

There remains to investigate the localization properties of \z >.The proba- 
bility density | < x\z > | 2 is maximal at a point Q determined by the equation 

dS 

y / 2hLu r Rez r (24) 



dQ 



4 



Eq.(24) gives an interpretation of Rez which coincides with the one of eq.(22) 
in the leading order of ?i.The probability density | < x\Ut\z > \ 2 is concentrated 
on the solution of the equation 

dS 



dQ r 



\J 2hcu r Re(exp(—iuj r t)z r ) (25) 



It can be seen from the definition (14) of Q that Q ~ x and S(Q) <~ Q 2 for 
the anharmonic oscillators in the weak coupling regime. So, in this sense the 
probability density | < x\U t \z > | 2 is concentrated on a deformed torus filled by 
the trajectories of the integrable system. 

Let us consider two examples. The first concerns a onedimensional non-negative 
potential V normalized in such a way that E r j = in cq.(ll). Then 



S cl (x) = V2ni Vv¥)dy (26) 



is uniquely defined (no problem with the square root ). Hence 



Q = exp{u^ I 1 dy) (27) 



To be more specific let 



V(y)^^-y 2 +gy i (28) 



We calculate the integral (26) 

S cl (x)=mLuJ dyy^l + cy 2 = ^(1 + cx 2 )^ (29) 
where c = . Now, 

Q = exp(J dy(y^l + cy 2 )- 1 ) = 2x{\ + + ex 2 )' 1 (30) 

As a second example let us consider a non-negative twodimcnsional rota- 
tionally symmetric potential V(r).Then, the classical Hamilton- Jacobi function 
W reads (in the plane coordinates (r, </>)) 

W(r,<j>)=h<j> + J dp(2mE cl (h,I 2 )-2mV(p)-I 2 p- 2 )^ (31) 

A solution of the "imaginary" Hamilton- Jacobi equation (11) is obtained in the 
limit E c i(Ii 7 I 2 ) — » 0. For a non- negative potential this limit means I\ — > and 
h — > 0. It follows that 

S cl (r,^)= J r dp(2mV(p))^ (32) 
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and 

Qi = exp(i(f> + ujiu(r)) (33) 
Q 2 = exp(2uj 2 u(r)) 

where 

u(r) = J dp(2mV{p))-? (34) 

and tOk are denned in eq.(19). 

So far we have restricted ourselves to the limit h — > . The corrections to 



the formulae (8)-(10) and (13) are discussed in our papers [JlT | (12 1 .In particular, 
in ref. jij) we bound the growth in time of the 0(K) term in eq.(8) by \^hi. 
We specify also some conditions which ensure that the 0{K) term is bounded 
uniformly in time. On the basis of these results we estimate that the wave 
packet moves along the trajectory (25) and does not change its shape at least 
for a time of order i . In general, it is known (|] (see also a rigorous formulation 
in ref.(8|) that a particle in a state belonging to the subspace of the Hilbert space 
corresponding to the discrete spectrum remains localized in a bounded region 
for an arbitrarily large time. Moreover, from the quasiperiodicity it follows 
that the wave packet recovers its shape after a sufficiently large time (see an 
estimate on this time in ref. ||). Our claim is that the coherent states defined in 
this Letter do not lose their shapes during a time interval large in comparison 
to the classical periods . 
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